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Abstract It has long been speculated that biological evolution was influenced by ultraviolet radiation 
(UVR) reaching the Earth's surface, despite imprecise knowledge of the timing of both UVR flux and 
evolutionary events. The past strength of Earth's dipole field provides a proxy for UVR flux because of its role 
in maintaining stratospheric ozone. The timing of Quaternary evolutionary events has become better 
constrained by fossil finds, improved radiometric dating, use of dung fungi as proxies for herbivore 
populations, and improved ages for nodes in human phylogeny from human mitochondrial DNA and Y 
chromosomes. The demise of Neanderthals at ~41 ka can now be closely tied to the intensity minimum 
associated with the Laschamp magnetic excursion, and the survival of anatomically modern humans can be 
attributed to differences in the aryl hydrocarbon receptor that has a key role in the evolutionary response 
to UVR flux. Fossil occurrences and dung-fungal proxies in Australia indicate that episodes of Late 
Quaternary extinction of mammalian megafauna occurred close to the Laschamp and Blake magnetic 
excursions. Fossil and dung fungal evidence for the age of the Late Quaternary extinction in North America 
(and Europe) coincide with a prominent decline in geomagnetic field intensity at ~13 ka. Over the last ~200 
kyr, phylogeny based on mitochondrial DNA and Y chromosomes in modern humans yields nodes and 
bifurcations in evolution corresponding to geomagnetic intensity minima, which supports the proposition 
that UVR reaching Earth's surface influenced mammalian evolution with the loci of extinction controlled by 
the geometry of stratospheric ozone depletion. 


Plain Language Summary The strength of Earth's magnetic field in the past, recorded by rocks 
and sediments, provides a proxy for past flux of ultraviolet radiation (UVR) to Earth's surface due to the role 
of the field in modulating stratigraphic ozone. About 40,000 years ago, mammalian fossils in Australia 
and Eurasia record an important die-off of large mammals that included Neanderthals in Europe. In the 
Americas and Europe, a large mammalian die-off appears to have occurred ~13,000 years ago. Both die-offs 
can be linked to minima in Earth's magnetic field strength implying that UVR flux variations to Earth's 
surface influenced mammalian evolution. For the last ~200,000 years, estimates of the timing of branching 
episodes in the human evolutionary tree, from modern and fossil DNA and Y chromosomes, can be linked to 
minima in field strength, which implies a long-term role for UVR in human evolution. New fossil finds, 
improved fossil dating, knowledge of the past strength of Earth's magnetic field, and refinements in the 
human evolutionary tree, are sharpening the focus on a possible link between UVR arriving at the Earth's 
surface, magnetic field strength, and events in mammalian evolution. 


1. Introduction 


The apparent spacing of mass extinction events at long (~26 Ma) repeat times (e.g., Raup & Sepkoski, 1986), 
and the supposed role of geomagnetic polarity reversal in extinction (e.g., Raup, 1985), has received 
intermittent attention over the last 50 years, since early studies of Quaternary radiolarian evolution and 
polarity reversal in deep-sea sediments (Hays, 1971). The reader is referred to Aubry et al. (2009) for 
abbreviations denoting geological time in the past (Ma for millions of years ago and ka for thousands of years 
ago) and equivalent durations (Myr and kyr, respectively). Linkage between polarity reversal (with its 
concomitant low field intensity) and extinction or speciation has not gained significant traction, perhaps 
because of uncertainties in the polarity timescale itself, and in the chronology of extinction/speciation 
outside of the few well-documented mass extinctions. On the other hand, we now know from Quaternary 
studies that although polarity reversals coincided with relative paleointensity (RPI) minima, intervals 
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between polarity reversals are also characterized by numerous RPI minima, some of which coincide with 
magnetic excursions (see Laj & Channell, 2015). The chronology of both the RPI record and the paleontolo- 
gical record remains poorly constrained, even for the Quaternary, such that a linkage between extinction 
and RPI minima cannot be ruled out. 


The geomagnetic field helps to preserve stratospheric ozone, as well as atmospheric composition, density, 
and oxygen levels that are vital to Earth's biosphere (Wei et al., 2014). The field shields Earth from galactic 
cosmic rays (GCRs) and solar wind, and from harmful ultraviolet radiation (UVR) that affect the function 
of living systems (Belisheva et al., 2012; Mendoza & de La Pefia, 2010). The demise of the Martian magnetic 
field, several billion years ago, is widely believed to have been the root cause for the near disappearance of 
the Martian atmosphere and the resulting dramatic change in the Martian environment from one featuring 
surface water and aqueous sedimentation to its present relative inactivity and sterility. The explosion of life 
in the Early Cambrian period at ~530 Ma has been associated with growth of Earth's inner core, the sup- 
posed strengthening of the dipole geomagnetic field, and the resulting thickening of Earth's atmosphere 
(Doglioni et al., 2016), although there is little evidence for strengthening of the geomagnetic field at this 
time (e.g., Biggin et al., 2015). On the other hand, the Late Ediacaran and Early Cambrian periods (~550 
and ~530 Ma, respectively) may have been times of unusually high polarity reversal frequency 
(Bazhenov et al., 2016; Pavlov & Gallet, 2001), although precise estimates of reversal frequency are elusive 
due to poorly constrained age control in stratigraphic sections where the reversals were recorded. Meert 
et al. (2016) proposed that high reversal frequency (up to ~20 reversals/Myr) at this time would have been 
associated with low geomagnetic field intensity that therefore lowered shielding from UVR, which created 
an evolutionary advantage for burrowing and shelled organisms. These proposals for the role of the geo- 
magnetic field in evolution are controversial partly because of poor knowledge of the state of the geomag- 
netic field 500-550 million years ago. Oxygenation of the oceans and atmosphere after the Gaskiers 
glaciation at ~580 Ma (Canfield et al., 2007) may have been the principal driver of the Early Cambrian 
explosion of life, both through oxygen levels at Earth's surface and increased UVR shielding through 
enhanced stratospheric ozone concentrations. 


Several strategies in modern organisms reflect the evolutionary impact of UVR. Behavioral adaptations to 
UVR include vertical water-column migration in aquatic organisms, the presence of UVR-screening pig- 
mentation decreasing with water depth, and complete disappearance of pigments for deep-water and 
cave-dwelling animals (e.g., Hessen, 2008). The red coloration of alpine plankton and the “red sweat” of 
the hippopotamus (Saikawa et al., 2014) are examples of evolutionary adaptation to high UVR at altitude 
and at low latitudes, respectively. UVR causes two classes of DNA lesions: cyclobutane pyrimidine dimers 
and 6-4 photoproducts. Both lesions distort DNA structure, introducing bends or kinks and thereby 
impeding transcription and replication (e.g., Branze & Foiani, 2008; Clancy, 2008). Relatively flexible areas 
of the DNA double helix are most susceptible to damage. One “hot spot” for UV-induced damage is found 
within a commonly mutated oncogene TP53 (Benjamin & Ananthaswamy, 2006), which in normal func- 
tion has an important role in tumor suppression. At low concentrations, reactive oxygen species (ROS) play 
vital roles during mutagenic activity in response to pathogen attack. Higher concentrations of ROS pro- 
duced by UVR give rise to oxidative stress where ROS attack DNA bases and the deoxyribosyl backbone 
of DNA (see MacDavid & Aebisher, 2014). Production of antioxidant enzymes neutralizes ROS, and ROS 
modulation is controlled by the aryl hydrocarbon receptor (AhR) that plays a key role in 
mammalian evolution. 


The consequences of ionizing radiation associated with GCRs and solar particle events (solar wind) for 
human health have received attention in recent years in an effort to evaluate the health effects of future 
space travel outside Earth's protective magnetosphere (e.g., Delp et al., 2016). Earth's atmosphere is opaque 
to all but the highest energy GCRs, and, together with the geomagnetic field, serves to shield Earth's surface 
from GCRs. The intensity of UVR arriving at Earth's surface decreases with increasing latitude and is atte- 
nuated by stratospheric ozone (O3) that acts as a sink for UVR. The geomagnetic field plays an important role 
in preserving the atmosphere, including stratospheric ozone that would otherwise be stripped away by solar 
wind and GCRs (Wei et al., 2014). UVR triggers dissociation of oxygen molecules (O2) into oxygen radicals 
that combine to form stratospheric ozone that in turn absorbs UVR as it splits into oxygen atoms. Certain 
ozone-depleting agents (such as nitrogen oxides) are produced naturally by energetic particle precipitation 
from solar wind, particularly during solar proton events, and therefore, times of low geomagnetic field 
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strength lead to higher ozone depletion (Randall et al., 2005, 2007). Atmospheric modeling implies substan- 
tial increases in hydrogen and nitrogen oxide concentrations due to enhanced ionization by GCRs during the 
Laschamp excursion, with significant decrease in stratospheric ozone particularly at high latitudes (Suter 
et al., 2014). Modeling of ozone depletion during polarity reversals, based on a geomagnetic field intensity 
~10% of the present value, leads to enhanced UVR flux at the Earth's surface, particularly at higher latitudes, 
that is 3-5 times that resulting from the anthropogenic ozone hole (Glassmeier & Vogt, 2010; Winkler et al., 
2008). Prior to anthropogenic emission of ozone-depleting chlorofluorocarbons and halons, energetic parti- 
cle precipitation at times of low geomagnetic field strength played an important role in ozone depletion. A 
well-defined nitrate peak, together with a broader 10Be peak, are associated with low field strength at the 
time of the Laschamp magnetic excursion (~41 ka) in the EPICA-Dome C Antarctic ice core (Traversi 
et al., 2016), which indicates that geomagnetic shielding played a role in the production of both cosmogenic 
isotopes (such as *°Be) and ozone-depleting nitrogen compounds. It is noteworthy that bacterial UVR 
proxies in sediments from Lake Reid (Antarctica) imply more than three times higher UVR flux during part 
of the last glacial than during the Holocene (Hodgson et al., 2005). UVR exposure affects the early stage of 
life in modern marine plankton, and plankton-benthos coupling in coastal waters (e.g., Hernandez 
Moresino et al., 2011). Furthermore, UVR plays a role in photosynthesis (e.g., Hollosy, 2002) and can cause 
changes in vegetation and habitat modification. 


The aim of this paper is to review the record of geomagnetic field intensity over the last ~300 kyr, and com- 
pare this record with the fossil record of extinction in terrestrial mammals (including Neanderthals), and 
with nodes in hominin phylogeny determined from mitochondrial DNA (mtDNA) and Y chromosomes. 


2. The Geomagnetic Field 


Knowledge of the Holocene geomagnetic field (i.e., over the last ~12 kyr) has been based on models built 
from archaeomagnetic, lava, and sediment data (Constable et al., 2016; Korte et al., 2011; Pavón-Carrasco 
et al., 2014). Beyond the Holocene, geomagnetic field strength during the Quaternary, over the last ~2.5 
Myr, has been acquired primarily from RPI data from marine sediments recovered by deep-sea drilling 
(e.g., Channell et al., 2009, 2018; Laj et al., 2000, 2004; Valet et al., 2005; Xuan et al., 2016; Ziegler et al., 
2011). Absolute paleointensity data from lavas are hampered by unknown time gaps between lava flows 
and inadequate age control in young (Quaternary) lavas. RPI data from sediments have been acquired by 
normalizing the natural remanent magnetization (NRM) intensity by a laboratory-acquired magnetization 
designed to activate the same population of magnetic grains that carry the NRM, thereby compensating 
for variations in concentration of NRM-carrying grains throughout the sedimentary section. The 
laboratory-applied normalizers are typically anhysteretic remanent magnetization (ARM) and/or isother- 
mal remanent magnetization (IRM). ARM is acquired in a decreasing alternating field (AF, with a peak 
AF of ~100 mT) with a weak direct current bias field (typically 50 uT), and IRM is acquired in a strong direct 
current field (up to ~1 T). The appropriate normalizer is usually chosen such that its coercivity (response to 
AF demagnetization) closely matches that of the NRM. A typical RPI proxy comprises the slope of NRM 
intensity versus ARM (and/or IRM) intensity determined over a particular peak demagnetization field range 
(such as 20-80 mT). The definition of the slope, often determined at 1-cm spacing down-core, can be gauged 
by determining its linear correlation coefficient (e.g., Xuan & Channell, 2009). 


There is substantial agreement among Late Quaternary stacks of global sedimentary RPI data (Channell 
et al., 2009; Valet et al., 2005; Xuan et al., 2016; Ziegler et al., 2011) that are independently supported by 
1Be/?Be data from marine sediments (e.g., Simon et al., 2016) due to the role of the geomagnetic field in 
modulating cosmogenic isotope (‘°Be) production. Agreement is, however, poor among stacks of sedimen- 
tary RPI data and /°Be/*Be data in the 0-40 ka interval (see Figure 1 of Channell et al., 2018) that is attributed 
to subtle and often unrecognized drilling disturbance in poorly consolidated uppermost sediments recovered 
from the ocean floor. A recently published RPI stack for 10-40 ka based on high-deposition-rate sediments 
from the Iberian Margin and elsewhere (Figure 1) is supported by paleointensity estimates from revised cal- 
culation of "Be flux in Greenland ice cores (Channell et al., 2018) using the GICCOS Greenland ice-core age 
model (Svensson et al., 2008). Models and stacks covering the same 0- to 40-ka interval that use RPI data 
from lower sedimentation rate sequences (e.g., Panovska et al., 2018) cannot resolve the detail that is 
revealed by higher sedimentation rate sequences and by Be flux in ice cores. 
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Figure 1. Relative paleointensity (RPI) calibrated to virtual axial dipole moment (VADM) for the interval 0-45 ka (a) and 
40-300 ka (b). (a) Overall RPI stack of 23 sedimentary records (red) with standard error (20) compared with the 10Be-based 
VADM from Greenland ice cores (blue) filtered with cutoffs of 1/3,000 years (dark blue) and 1/500 years (dashed light 
blue; Channell et al., 2018). Also shown are Holocene VADM models of Constable et al. (2016; green) and Pavón-Carrasco 
et al. (2014; black). (b) PISO RPI stack for 40-300 ka with standard error (Channell et al., 2009). The yellow shading 
indicates RPI minima in the stack that have been associated with directional magnetic excursions. See text for references 
that document the labeled excursions. 


When adequately recorded, “magnetic excursions” are manifested in both lavas and sediments as antipodal 
magnetization directions defining short-lived polarity reversal, and they occupy minima in RPI records (Laj 
& Channell, 2015). At least seven magnetic excursions have been documented in lavas and sediments of the 
last 300 kyr (Figure 1) with ~13 in the Brunhes Chron (last 775 kyr). Directional magnetic excursions have 
been named after the location where they were initially recorded such as Mono Lake (34 ka), Laschamp (41 
ka), Blake (~115 ka), Iceland Basin (190 ka), and Pringle Falls (211 and/or 238 ka). Although the first record 
of a magnetic excursion (the Laschamp excursion) was published over 50 years ago (Bonhommet & Babkine, 
1967), magnetic excursions remain controversial mainly because the aberrant magnetization directions that 
define them have submillennial to millennial duration and are therefore fortuitously recorded in geological 
archives. Although the Laschamp excursion has now been recorded in scores of globally distributed sedi- 
mentary sequences (e.g., Channell et al., 2017; Collins et al., 2012; Evans et al., 2007; Laj et al., 2000, 2006; 
Lund et al., 2005; Mazaud et al., 2002), the Laschamp records represent a small fraction of the total number 
of paleomagnetically studied sedimentary sequences that cover the Laschamp interval. The submillennial 
duration of the excursion means that recordings of the excursion are generally restricted to sequences with 
mean sedimentation rates >10-15 cm/kyr that are unusual in deep-sea sediments (see Roberts & 
Winklhofer, 2004). Lava records of the Laschamp are restricted to the Massif Central in France (e.g., Laj 
et al., 2014), the Auckland volcanic field (Cassata et al., 2008), and New Zealand's Mt. Ruapehu (Ingham 
et al., 2017), that just happened to have erupted with high frequency during this brief period. 
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The age of the Laschamp excursion, and the low geomagnetic field strength associated with it (Figure 1), is 
based on correlation of deep-sea cores that record the Laschamp excursion to ice-core chronologies (e.g., 
Channell et al., 2017; Laj et al., 2000; Nowaczyk et al., 2012), by 4 Ar/PPAr age determinations in lavas from 
the Massif Central (e.g., Laj et al., 2014), by U/Th ages in a speleothem from Missouri (Lascu et al., 2016), and 
by the age of the related cosmogenic-isotope flux peak in ice cores (Baumgartner et al., 1998; Svensson et al., 
2006, 2008; Traversi et al., 2016; Wagner et al., 2000; Yiou et al., 1997). Current estimates of the duration and 
midpoint age of the Laschamp directional excursion are <1 kyr and 41 ka, respectively. A less abrupt geo- 
magnetic field intensity decrease (over several kyr) defines a paleointensity minimum that brackets the 
directional excursion. 


Directional records of the Iceland Basin excursion at ~190 ka have been recovered from numerous deep-sea 
sediment cores from the Northern Hemisphere (see review in Channell, 2014), in cores from Lake Baikal 
(Oda et al., 2002), and two sites from the South Atlantic ocean (Channell et al., 2017; Stoner et al., 2003). 
The age of the excursion has been established by correlation to marine oxygen isotope records, and to ice 
core records. There are no records of the Iceland Basin excursion from well-dated lava sequences. 


Apart from the Laschamp and Iceland Basin excursions, magnetic excursions during the last 300 kyr are con- 
troversial because of the paucity of records and/or poor age control in available records. On the other hand, 
RPI records and hence recordings of RPI minima associated with the excursions are numerous, presumably 
because RPI minima associated with excursions are manifest over longer timescales than the associated 
directional excursions. 


The Mono Lake excursion at ~34 ka has been recorded in both deep-sea sediments and lava flows (Cassata 
et al., 2008; Channell, 2006; Kissel et al., 2011; Laj et al., 2014; Negrini et al., 2014) as well as at Mono Lake in 
California (e.g., Benson et al., 2003), although the age at the Mono Lake type-location is controversial (e.g., 
Cassata et al., 2010; Kent et al., 2002). The age of a RPI minimum associated with the Mono Lake excursion 
can be estimated from a cosmogenic-isotope flux maximum in Greenland ice cores (Muscheler et al., 2005; 
Wagner et al., 2000). 


Although the older of the two Blake excursions (100 and 120 ka) was first recorded almost 60 years ago 
(Smith & Foster, 1969), there have only been a handful of observations of either Blake excursion since 
(see Laj & Channell, 2015). The older Blake excursion has, however, been observed in diverse media: marine 
sediments (e.g., Tric et al., 1991), Chinese loess (e.g., Zhu et al., 1994), and a speleothem from Spain (Osete 
et al., 2012). The younger “Blake” excursion, also known as the Skalamaelifell excursion, has been observed 
in Icelandic and Italian lavas (Jicha et al., 2011; Singer et al., 2014). The age of the Pringle Falls excursion (at 
211 ka and/or 238 ka) is not yet settled due to a lack of records in well-dated sequences (Herrero-Bervera 
et al., 1994; Laj & Channell, 2015; Singer, 2014; Singer et al., 2008). 


An apparent RPI minimum at ~13 ka appears in some individual RPI records, as a notch in RPI stacks 
(Figure 1), and in virtual axial dipole moment (VADM) proxies from 10Be flux in Greenland ice cores 
(Channell et al., 2018). In one core from the Iberian Margin (MD01-2444), a pronounced RPI minimum at 
~13 ka is associated with a directional magnetic excursion (Channell et al., 2013). The apparent RPI mini- 
mum at ~13 ka in core MDO1-2444 coincides with the old-end of the 0-14 ka Holocene paleointensity model 
of Pavon-Carrasco et al. (2014). The RPI stack of Channell et al. (2018) contains a prominent field intensity 
peak (VADM ~13 x 10° Am”) around 18-17 ka followed by a field collapse between 16 and 14 ka culminat- 
ing ina minimum VADM value (~7 x 10?” Am?) that has been smoothed by stacking (Figure 1). The VADM 
decrease at 17-13 ka corresponds to a rate of decrease (12.1 nT/year) that is comparable to the rate (12.7 
nT/year) observed in the 1Be-derived field at the same time (Channell et al., 2018) and to the rates observed 
during the Laschamp/Mono Lake excursions (Laj & Kissel, 2015). 


To provide a reference template for geomagnetic field intensity over the last ~300 kyr, we combine the RPI 
stack for 14-45 ka (Channell et al., 2018) with the Holocene model for 0-14 ka (Pavén-Carrasco et al., 2014) 
and with the PISO paleointensity stack (Channell et al., 2009) beyond 45 ka. This paleointensity template 
was calibrated to VADM by aligning it with Holocene archaeomagnetic models (Constable et al., 2016; 
Korte et al., 2011; Pavén-Carrasco et al., 2014) and assuming a value of ~1.5 x 10°? Am? for the VADM mini- 
mum at the Laschamp excursion (Laj et al., 2014). This VADM template provides a proxy for geomagnetic 
intensity variations that is compared with events in mammalian evolution, in order to test potential linkages 
between geomagnetic field intensity variations and Late Quaternary mammalian evolution. 
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3. Late Quaternary Extinctions 


Causes of extinction of mammalian megafauna (adult weight >45 kg) during the Late Quaternary, the so-called 
Late Quaternary Extinction (LQE), have been the subject of prolonged debate (e.g., Koch & Barnosky, 2006; 
Martin, 1967; Stuart, 2015). Prior to ~13 ka, the mammal assemblage of the Americas included large-bodied ani- 
mals such as mammoths, horses, camels, saber-tooth cats, and the short-faced bear. Extinction was total for 
mammals larger than 1,000 kg, >50% for size classes between 32 and 1,000 kg, and ~20% for those between 
10 and 32 kg (Koch & Barnosky, 2006). Within a short time window, >150 species were lost in the Americas, 
including all mammals over ~600 kg. An analogous catastrophic size-controlled LQE affected 14 of 16 
Australian mammalian genera; however, extinction ages are at least ~30 kyr older than in North America. 
Although fewer species were affected by the LQE in Africa and Eurasia, a similar size-biased extinction has 
been observed with end-Pleistocene (~13 ka) ages being prominent. Current explanations for the LQE in 
North America and Australia involve a combination of two hypotheses: climate change, and “overkill” by 
human hunting, modulated by the knock-on effect of herbivore extinction on the environment and on the sur- 
vivability of other groups (e.g., Owen-Smith, 1987). Although “overkill” was originally used to explain North 
American extinctions (Martin, 1967), a forerunner of the hypothesis was popular in 19th century Europe, before 
it was eventually abandoned as archaeological evidence for human migration showed little evidence for the 
impact of human hunting on the LQE. Grayson and Meltzer (2003) argued that island settings (e.g., New 
Zealand or the West Indies), where human hunting and habitat degradation can be unequivocally associated 
with extinction, should not be the model for continental extinctions. Extraterrestrial impact as a contributing 
cause for the LQE in North America, and for the onset of the Younger Dryas cold period (Firestone et al., 
2007), has not been supported by subsequent analyses (Holliday et al., 2014; Pinter et al., 2011). 


In North America, the brief (~200 year) duration of Clovis-tool finds at ~13 ka (e.g., Waters & Stafford, 2007) 
is usually associated with rapid dispersal of modern humans across North America and is closely contem- 
poraneous with the LQE peak. On the other hand, the lack of mammalian kill sites in the Clovis record 
argues against a direct linkage between Clovis technology and “overkill.” Humans were present in North 
America at least several kiloyear prior to the Clovis horizon (Gilbert et al., 2008; Waters et al., 2011), and per- 
haps prior to ~24 ka (Bourgeon et al., 2017). 


According to Faith and Surovell (2009), the LQE in North America was abrupt and requires a mechanism 
capable of wiping out ~35 genera across the continent in a “geological instant” in the 13.8- to 11.4-ka interval 
(Figure 2), with the spread in last appearances being largely explained by an incomplete fossil record and the 
resulting Signor-Lipps effect (Signor & Lipps, 1982). Abrupt versus staggered megafaunal extinction at the 
LQE is central to determination of cause. In their “continental simulation,” Faith and Surovell (2009) deter- 
mined the empirical probability (3.4%) of observing a terminal Pleistocene (10-12 ka) age from 1955 strati- 
graphic occurrences (from 31 genera) of which 66 taxa (from 16 genera) yield terminal Pleistocene ages, 
assuming that all occurrences are equally likely to receive a terminal Pleistocene age. The simulation ran- 
domly assigned pre-12-ka or post-12-ka ages to all 1955 observations based on this probability (3.4%). The 
total number of genera that received a terminal Pleistocene age was tallied for each of 10,000 simulations 
to determine the probability of observing 16 or fewer terminal Pleistocene genera. The authors concluded 
that the observed pattern is consistent with synchronous (i.e., 10-12 ka) extinction for all 31 genera. 


Bradshaw et al. (2012) proposed a Gaussian-resampled, inverse-weighted McInerny (GRIWM) approach, 
which weights observations inversely according to their temporal distance from the last observation of a con- 
firmed species occurrence, and samples radiometric ages from the underlying probability distribution. In 
Figure 2, we show GRIWM estimates of continent-wide European extinctions from the fossil record aided 
by DNA analyses (Cooper et al., 2015), excluding regional disappearances. An extinction age estimate in 
North America for Arctodus simus (the short-faced bear) at 10.8 ka (Schubert, 2010) and the onset of popula- 
tion decline of Bison priscus (the steppe bison) at ~37 ka (Shapiro et al., 2004) are included in Figure 2. Note 
that the horse and woolly mammoth (Mammuthus primigenius) persisted in interior Alaska until ~10.5 ka 
(Haile et al., 2009), and the woolly mammoth survived on St. Paul Island (Alaska) until ~5.6 ka (Graham 
et al., 2016). Zazula et al. (2014) pointed out that the American mastodon (Mammut americanum) occupied 
eastern Beringia (Alaska/Yukon) during the last interglacial before its range contracted southward at the 
onset of glacial conditions at ~75 ka. The range of the species appears to have expanded northward again 
as interglacial conditions returned at the end of the Pleistocene, before extinction of the species at ~11.5 
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Figure 2. Geomagnetic field intensity (VADM) for the last 70 kyr compared with timing of continent-wide extinction events for genera from Australia (Saltré et al., 
2016), Eurasia (Cooper et al., 2015), and North America (Faith & Surovell, 2009), including extinction of Arctodus simus (Schubert, 2010) and population decline for 
Bison priscus (Shapiro et al., 2004). The solid and dashed bars in the Eurasia record refer to fossil and DNA derived events, respectively. Last appearance 

datums for extinct North American genera (Faith & Surovell, 2009) were recalibrated using Calib 7.1 (Stuiver et al., 2018) with 20 errors for entries that pass 
selection criteria established by Meltzer and Mead (1985). Dung fungi (Sporormiella and Podospora) in sedimentary archives, proxies for herbivore population, are 
for North America (Gill et al., 2009) and Australia (Johnson et al., 2015; van der Kaars et al., 2017). Geomagnetic field intensity during the last 70 kyr (red) 
combines VADM for the last 14 kyr (Pavon-Carrasco et al., 2014), VADM for 14-45 ka (Channell et al., 2018), and the PISO VADM record (Channell et al., 2009) for 
45-70 ka, The light blue band represents the VADM estimate from *°Cl-flux in Greenland ice-cores (Muscheler et al., 2005). The dark blue line corresponds to the 


ice-core 


10Be-derived VADM record (Channel « et al., 2018). The background shades of blue are scaled to the maximum (darker blue color) and minimum 


(white) VADM intensity in the 1.35-13.27 x 10°? Am? range. 


ka (10,000 "ÍC years before present [BP]). Zazula et al. (2014) posed the question: why was this species 
stopped in its tracks when favorable conditions beckoned in Beringia? 


An important proxy for herbivore population is the abundance in sedimentary sequences of coprophilous 
(dung) fungal spores, such as Sporormiella. The proxy was first proposed over 30 years ago (Davis, 1987), 
requires careful interpretation and laboratory techniques (e.g., van Asperen et al., 2016), but provides a mea- 
sure of herbivore population independent of the bone-fossil record. Lake sediments in New York and 
Indiana imply a decline in Sporormiella spores beginning at 14.8 ka (Figure 2) that falls below the 2% thresh- 
old by 13.7 ka (Gill et al., 2009). This result has been closely replicated at Silver Lake (Ohio) where 
Sporormiella decline was dated at 13.9 ka (Gill et al., 2012). Importantly, the Sporormiella decline at these 
sites predates Younger Dryas cooling, and concurrent changes in the pollen record, and immediately pre- 
dates a marked charcoal deposition increase, implying that herbivore decline and resulting landscape 
changes provide an explanation for subsequent (natural) landscape burning. The onset of the demise of 
North American herbivores at ~14.5 ka (Gill et al., 2009) lies within the Bolling-Allerad warm period with 
the Younger Dryas cold period beginning ~2 kyr later (e.g., Deplazes et al., 2013). 


The South American LQE was even more profound than that in North America (Barnosky & Lindsey, 2010; 
Koch & Barnosky, 2006), with the loss of 50 megafaunal genera (~83%). Robust dates are scarce for the South 
American LQE, although it appears that many taxa were lost near the Pleistocene-Holocene boundary 
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(Barnosky & Lindsey, 2010). Sporormiella decline in lake sediments from SE Brazil implies herbivore popu- 
lation collapse at ~12 ka (Raczka et al., 2018). 


In Northern Eurasia, nine genera (35%) were lost during the LQE. Available age data are consistent with a 
two-phase extinction in the 45- to 35-ka and 15- to 10-ka intervals (Koch & Barnosky, 2006). Up to 50% of 
worldwide megafaunal extinctions at 15-10 ka apparently occurred in Northern Eurasia (Cooper et al., 
2015), but the extinction pattern is more complex than in North America with megafaunal range contrac- 
tions culminating in extinction for some species but not others (Stuart, 2015). In Figure 2, we plot 
continent-wide Eurasian megafaunal extinction events from Cooper et al. (2015). Several well-studied spe- 
cies disappeared continent-wide at ~26-31 ka (Figure 2); hence, their last appearances significantly postdate 
the Laschamp excursion (at ~41 ka). Post-Laschamp extinction for Crocuta crocuta (spotted hyaena) and 
Crocuta spelaea (cave hyaena) at ~26 ka were, however, preceded by severe range contraction from Asia into 
Europe (Stuart & Lister, 2014). Fossils of Ursus spelaeus (cave bear) also indicate E to W range contraction 
before extinction at ~26 ka with abrupt population decline, based on DNA analyses, after 50 ka (Baca 
et al., 2016; Stiller et al., 2010; Stiller et al., 2014). 


In Africa, at least 24 species and ~10 genera of mammals became extinct in the 13- to 6-ka interval, represent- 
ing 25% of Pleistocene African megafauna (Faith, 2014). Species-level extinction was, again, most intense for 
larger-bodied megafauna (Koch & Barnosky, 2006). The African LQE was considered to have been less 
severe than elsewhere, accounting for the relatively rich diversity of modern African megafauna. On the 
other hand, the number of extinct African species dated to the last 100 kyr exceeds the number in Europe 
and matches the number in Australia, and is only surpassed by the LQE in the Americas (Faith, 2014). In 
east Africa alone, the number of securely dated latest Pleistocene mammal extinctions has risen from two 
to seven in the last decade, with most being of grazers associated with open habitats (Faith, 2014). 


Estimation of extinction ages for Australian megafauna (and for some Eurasian genera) is complicated by 
the majority of last appearances being at or beyond the practical range of radiocarbon dating (i.e., >40 
ka). The LQE in Australia was apparently catastrophic for large mammalian megafauna, with the complete 
loss of all animals heavier than ~100 kg. Fourteen of 16 genera of Pleistocene mammalian megafauna disap- 
peared, together with all megafaunal reptiles (6 genera), in the vicinity of, or prior to, ~40 ka (Figure 2). Ten 
Australian genera disappeared in the 44- to 35-ka interval based on a variety of frequentist statistical meth- 
ods (including GRIWM) to determine extinction ages for 16 megafaunal genera (Figure 2; Saltré et al., 2016). 
The mass extinction of megafauna at this time, including the largest-known (~3,000 kg) marsupial 
(Diprotodon), has been linked with climate variability and aridity (e.g., Wroe et al., 2013), although this link- 
age has been disputed (e.g., Saltré et al., 2016), often in favor of human predation or “overkill” (e.g., Brook & 
Johnson, 2006; Johnson et al., 2016; Miller et al., 2016; van der Kaars et al., 2017). It is noteworthy that the 
extinction age for the ~200-kg flightless bird Genyornis newtoni at ~35 ka (Figure 2; Saltré et al., 2016) is 
younger than the ~43-ka estimate given by Miller et al. (2016) based on dated eggshell fragments. Even if 
final extinction was delayed until ~35 ka, the population of Genyornis newtoni crashed close to the time of 
the Laschamp excursion (~41 ka) based on the egg-shell data (Miller et al., 2016), although some egg-shells 
attributed to Genyornis newtoni may be from other species (Grellet-Tinner et al., 2016). 


At Lynch's Crater (NE Queensland), an abrupt decline in dung fungi including Sporormiella (Figures 2 and 
3) implies abrupt demise of large Australian herbivores at 40-44 ka Johnson et al., 2015). An abrupt increase 
in charcoal lags Sporormiella decline by ~100 years, and evidence for grasses and sclerophyll vegetation lags 
Sporormiella decline by ~300-400 years (Johnson et al., 2015; Rule et al., 2012). The charcoal-rich levels can 
be explained by natural lightning-induced biomass burning as a result of fuel buildup triggered by herbivore 
extinction (Johnson et al., 2016; Rule et al., 2012). Off the southern coast of Western Australia, marine core 
MD03-2614G records a sharp decline in Sporormiella in the 45- to 43-ka interval, relative to values recorded 
back to 140 ka (Figures 2 and 3; van der Kaars et al., 2017). 


A role for humans in the extinction of large animals in Australia remains popular (e.g., Brook & Johnson, 
2006; Johnson et al., 2016; Miller et al., 2016; Turney et al., 2008; van der Kaars et al., 2017), although the 
arrival of humans in Australia (Sahul) may have predated the LQE at ~40 ka by ~25 kyr (Clarkson et al., 
2017) but the arrival date is not unequivocal (O'Connell et al., 2018). There is no evidence for a spike in 
the human population in Australia at the time of the most prominent extinction event at ~40 ka, when 
the entire Australian human population may not have exceeded a few tens of thousands (Williams, 2013). 
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Figure 3. Geomagnetic field intensity (VADM) for the 35- to 50-ka interval compared to probability density functions for 
the timing of Neanderthal disappearance at different sites (blue) and for the end of the Mousterian tool industry (black; 
Higham et al., 2014), with dung-fungal proxies for large herbivore extinctions from NE Queensland (green; Johnson et al., 
2015) and southern Western Australia (open circles joined by purple line; van der Kaars et al., 2017). The yellow band 
represents the GLOPIS paleointensity record converted to VADM (Laj et al., 2004, 2014), the purple curve is the sedi- 
mentary VADM record with 2c errors, and the red line is the VADM reconstructed from 10Be flux in Greenland ice-cores 
(Channell et al., 2018). The gray shading indicates the age range associated with Neanderthal demise, and the temporal 
overlap between the two human groups (Neanderthals and Cro Magnon). The end of the Mousterian flint-tool industry 
has been dated at 39.260-41.030 ka (Higham et al., 2014), and 39.528-41.013 ka (Cooper et al., 2015), corresponding closely 
with low field intensity associated with the Laschamp geomagnetic excursion (39.7-41.9 ka; Laj et al., 2014). 


Tasmania had a land bridge to the Australian mainland during the last glacial, becoming an island in the 
early Holocene. The extinction of megafauna in Tasmania at ~40 ka does not correspond to climate or 
environmental change and has been associated with the late arrival of humans in the region (Turney 
et al., 2008). More recent results place the Tasmanian extinction of Protemnodon anak and other 
megafauna at ~41 ka, predating human arrival on the island at ~39 ka and hence precluding human 
involvement in the extinctions (Cosgrove et al., 2010; Lima-Ribeiro & Diniz-Filho, 2014). Extant smaller 
(more accessible) prey, particularly the common wombat (Vombatus ursinus) and the red-necked wallaby 
(Macropus rufogriseus), characterize the early archaeological kill-sites on the island (Cosgrove et al., 2010). 


Apart from Australian extinctions concentrated close to the Laschamp excursion at ~40 ka, fossils from the 
King's Creek Catchment (SE Queensland) indicate additional concentrations of megafaunal last 
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appearances at ~83, ~107, and ~122 ka (Price et al., 2011; Wroe et al., 2013). The older two dates (107 and 122 
ka) correspond to magnetic field intensity minima associated with the Blake excursions (Figure 1b). 


4. Neanderthal Extinction 


The extinction of Neanderthal represents one of the great puzzles of human evolution. Neanderthal and ana- 
tomically modern humans (AMHs) cohabited Western Europe for ~2-5 kyr, prior to ~39 ka, supporting the 
contention that competition may have contributed to the demise of Neanderthal (Higham et al., 2014). Brief 
cold and dry conditions in Europe associated with Heinrich Stadial (HS) 4 were proposed as an additional 
likely stressor on Neanderthal (Sepulchre et al., 2007). Analyses of Campanian Ignimbrite (CI) cryptotephra 
from archaeological sites in Greece and elsewhere in Eastern Europe and Libya indicated that the CI erup- 
tion occurred early in a dry period associated with HS4, postdated the end of the Middle Paleolithic and the 
Mousterian tool industry, and hence postdated the demise of Neanderthal (Lowe et al., 2012). In Black Sea 
sediment cores, CI tephra overlie, and therefore postdate, the Laschamp excursion (Nowaczyk et al., 2012). 


The extinction of Neanderthal and the demise of the Mousterian tool industry (Figure 3) can be placed at 
41,030-39,260 calibrated years before present (41-39 ka) with 95.4% probability (Higham et al., 2014). 
Cooper et al. (2015) estimated the extinction of Neanderthal at 41,227 calibrated years BP with a standard 
deviation of 219 years, and 39,528-41,013 calibrated years BP using the GRIWM method (Figure 2). 
Mousterian ages outside this range have been recorded at several locations in southern Iberia including 
Gorham's Cave in Gibraltar (Finlayson et al., 2006; Tzedakis et al., 2007), but these ages should now be disre- 
garded according to Higham et al. (2014). Recent findings cast doubt on the existence of Neanderthal after ~39 
ka and lead to a closer correspondence of the demise of Neanderthals with the Laschamp magnetic excursion 
and the associated brief interval of very low geomagnetic field intensity centered at ~41 ka (Laj et al., 2014). It is 
important to note that the IntCal13 radiocarbon calibration (Reimer et al., 2013) may be offset to older ages by 
~1 kyr in the vicinity of the Laschamp excursion, relative to ice-core chronologies (Muscheler et al., 2014). 


Valet and Valladas (2010) proposed that low magnetic field strength in the Laschamp/Mono Lake excursion 
interval (40-33 ka) was an important factor in Neanderthal demise. Why AMHs were not similarly affected 
has remained an open question considering that the two populations shared habitats for 2,600-5,400 years 
(Higham et al., 2014) or >5 kyr (Lowe et al., 2012). There is no evidence for differences in skin pigmentation 
between European AMHs and Neanderthals, and at least a fraction of Neanderthals apparently had the same 
pale skin and/or red hair observed in some modern humans (Lalueza-Fox et al., 2007). Natural skin pigmen- 
tation in humans mitigates the harmful effects of UVR, but its advantage is offset by the importance of sun- 
light for vitamin D3 synthesis. The skin protection factor of “red ochre” (hematitic iron oxides) is 
traditionally utilized by some African tribes (Rifkin et al., 2015) and has been used since at least the last 
interglacial (~120 ka) based on ochre coatings on strung beads and residues on storage shells from Africa 
(Hodgskiss & Wadley, 2017), SE Spain (Hoffmann et al., 2018), and Levantine Mousterian sites (Bar-Yosef 
Mayer et al., 2009). The mystery of AMH survival at the time of Neanderthal demise may have been resolved 
by the discovery of differences in amino acid substitution in an intracellular chemosensor (the AhR) for 
AMHs and for Neanderthals and other primates (Hubbard et al., 2016). 


5. The Role of the AhR 


Defense mechanisms against UVR include the production of quenching agents and antioxidant enzymes that 
neutralize ROS produced by UVR. The ROS modulation is controlled by the AhR intracellular chemosensor 
that plays a key role in the evolutionary response to UVR. Experimental results indicate an adaptive response 
of mitochondria to varying ROS levels under a phenomenon called mitohormesis (Becker et al., 2016). 


The primary role of AhR is to regulate the transcription of genes mediating responses to the biochemical and 
toxic effects of dioxins, polyaromatic hydrocarbons, and related compounds (Abel & Haarmann-Stemmann, 
2010). AhR is expressed in all skin cells and can be generated by UVR through an endogenous ligand formed 
in situ from an amino acid called tryptophan (Esser et al., 2009). UVR and the more harmful UVB (wave- 
length 290-320 nm) are absorbed by free tryptophan in the cytosol of epidermal cells, and AhR plays a key 
role in translocating UVR stress response to the nucleus (Fritsche et al., 2007; Tigges et al., 2014; Wei 
et al., 1999). Exposure to UVR, particularly UVB, generates highly mutagenic DNA photoproducts. The pro- 
cess initiates apoptosis and involves damage to nuclear DNA accompanied by mitochondrial dysfunction 
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(Frauenstein et al., 2013). There is a general consensus that the AhR of modern humans is implicated in 
DNA repair (Dittmann et al., 2016; Schreck et al., 2009), tumor suppression (Fan et al., 2010; Yu et al., 
2017), epidermal barrier function (Noakes, 2015), skin tanning response, and melanocyte homeostasis 
(Jux et al., 2011; Luecke et al., 2010). Phylogenetic analysis suggested that the ability of vertebrate AhR to 
sense xenobiotics was acquired at a late stage of evolution, implying that the driving force for evolutionary 
conservation of AhR lies not only in its role in xenobiotic metabolism but also in normal cell development 
(Hahn et al., 2017; Hao & Whitelaw, 2013). 


Hubbard et al. (2016) showed that the AhR variant in modern humans contains Val381 residue in the ligand- 
binding domain, while the AhR of Neanderthals, and a Denisovan individual, as well as nonhuman primates 
and other vertebrates (rodents) encode the ancestral Ala381 variant. The Val381 variant is fixed in the gen- 
ome of all modern humans as well as in the genome of the oldest (45 ka) AMH individual sequenced to date 
(Fu et al., 2014). Hubbard et al. (2016) suggested that the unique modification of AhR in AMHs led to signif- 
icant competitive advantage over their Neanderthal neighbors, due to decreased sensitivity in AMHs to tox- 
ins associated with fire-smoke, the effects of which may have been exacerbated by troglodytic lifestyles. 


Our focus here is on AhR involvement in the regulation of the skin responses to UVR, especially to harmful 
UVB radiation, and its modulation of the immune system (Agostinis et al., 2007; Esser et al., 2013; Rannug & 
Fritsche, 2006). UVB induces two signaling routes in mammalian cells: first, UVB is absorbed by nuclear 
DNA that results in generation of DNA photoproducts, and second, UVB activates cell-surface receptors 
(Merk et al., 2004). AhR plays an important role in skin integrity and immunity. AhR activation leads to 
transcriptional gene activation and is involved in the cutaneous stress response to UVR (Agostinis et al., 
2007; Dittmann et al., 2016; Navid et al., 2013; Schwarz, 2005) and alterations of gene expression (Dugo 
et al., 2012). Activation of AhR by UVB leads to signaling both to the nucleus and to cell membranes 
(Fritsche et al., 2007). The findings show that UVB irradiation affects cell surface receptors with subsequent 
activation of mitogen-activated protein kinases that in turn affect DNA in the nucleus. 


6. The Early Fossil Record of Homo sapiens 


From analyses of the fossil and stone-tool record over the last 250 kyr (Figure 4), Lahr (2016) proposed five 
transitions in the evolutionary history of Homo sapiens: (1) the origins of the species at 240-200 ka, (2) the 
first major expansion at 130-100 ka, (3) a period of dispersals at 70-50 ka, (4) a period of local/regional struc- 
turing of diversity at 45-25 ka, and (5) an early Holocene phase of significant extinction of hunter-gatherers 
and expansion of farmers (the Holocene Filter). 


Until recently, Member 1 of the (Omo) Kibish Formation of Ethiopia yielded the earliest known AMH cranial 
remains (Day, 1969). A volcanic tuff about 2 m below the level of the fossil finds has a 4 ar/? Ar weighted 
mean age of 196 + 2 ka (McDougall et al., 2005). An age of 196 ka for Member 1, combined with sedimento- 
logical evidence for rapid deposition, is consistent with Member 1 having been deposited synchronously with 
Mediterranean sapropel S7 (Brown et al., 2012; McDougall et al., 2008). Jebel Irhoud, Morocco, has been an 
important archaeological site since the 1960s when human fossils were found alongside Mousterian stone tools 
and were once dated at ~40 ka (Ennouchi, 1962). Recent fossil discoveries at this location support the presence 
of fossils with characteristics of AMHs (Hublin et al., 2017), and new luminescence age dating (Richter et al., 
2017) indicates that the fossils are considerably older, at ~300 ka, than the Ethiopian AMH finds at Omo Kibish 
(Figure 4). The finds in Morocco are now among the earliest known hominin fossils with AMH characteristics, 
in common with characteristics of the Florisbad fossil from South Africa discovered by T. F. Dreyer in 1932 and 
dated at ~260 ka (Griin et al., 1996). With such a sparse hominin fossil record, it is not possible to predict when 
AMHs first appeared, although the period from 300 ka to 200 ka in Africa appears to have been a critical time 
in development of AMH, and may constitute the evolutionary cradle of our species. 


7. Time to Most Recent Common Ancestor (TMRCA) From mtDNA and 
Y Chromosomes 


Thirty years ago, Cann et al. (1987) demonstrated that ethnically diverse surveys of modern mitochondrial 
DNA (mtDNA) are a major source of human evolutionary history. The advantages of mtDNA as an evolution- 
ary tool include faster mutation rates in mtDNA than in nuclear genes, and mtDNA is inherited maternally 
and does not recombine. Cann et al. (1987) demonstrated the African origin of the human mitochondrial 
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Figure 4. (a) Schematic evolutionary history of Homo sapiens modified after Lahr (2016). The fossil record includes key fossil specimens from Eliye Springs, 
Guomde, Florisbad, Herto, Omo Kibish, and hominins from Jebel Irhoud dated by Richter et al., (2017). The black/gray lines represent modern humans, the 
dark blue lines represent the Denisovans, and the light blue lines represent the Neanderthals. Neanderthal extinction has been drawn according to the timing of the 
end of the Mousterian tool industry at ~40 ka (Higham et al., 2014). (b) Estimates of the time to the most recent common ancestor (TMRCA) and phylogenetic 
nodes, from Y chromosomes (blue icons) and mitochondrial DNA (mtDNA; pink icons), according to different authors. For Scozzari et al. (2014), haplogroup B 
corresponds to the split between chromosomes found only in central-west Africa and chromosomes spread over sub-Saharan Africa, whereas the node at 115 ka 
marks the separation between African specific and all remaining haplogroups. Data from Wei et al. (2013; GENETREE-2 model) include the TMRCA of haplogroup 
A3 and three phylogenetic nodes considered by the authors of particular interest: DR (expansion of Y chromosomes following the out-of-Africa migration), FR 
(paleolithic male lineage expansion), and R1b (Neolithic modern European chromosomes). The data from Rito et al. (2013) refer to different nodes of the typical 
African haplogroup LO. (c) Geomagnetic VADM for the last 300 kyr (red band) as in Figure 2, using the PISO VADM record (Channell et al., 2009), to extend the 
record back to 300 ka. The light blue band and the dark blue line represent the VADM estimate C°cL and !°Be-derived, respectively) from ice-cores as in Figure 2. 
The yellow line represents the VADM determined from the GLOPIS paleointensity stack from marine sediments (Laj et al., 2004, 2014). Magnetic excursions 
corresponding to paleointensity minima are labeled (see Figure 1). The background shades of blue are scaled to the maximum (darker blue color) and minimum 
(white) VADM intensity in the 1.35-14.73 x 10°? Am? range. 


gene pool and estimated that mtDNA stems from an African mitochrondrial “Eve” who lived approximately 
200,000 years ago. Initially, this conclusion met with considerable resistance, however, increasing numbers 
of studies, including work on Y-chromosomes, imply that the Eve hypothesis is substantially correct. There 
is now broad consensus for the “out-of-Africa” hypothesis, whereby modern humans appeared at ~200 ka 
in Africa and spread throughout the continent before dispersing across the globe, although the exact chronol- 
ogy and nature of the population divergence remain unclear (Zhou & Teo, 2016). Dispersal has resulted in the 
occupation of a wide variety of habitats with selection in response to specific ecological pressures. A complete 
understanding of adaptation depends on a description of the genetic mechanisms and selective history that 
underlies heritable traits (e.g., Radwan & Babik, 2012), and the signatures of natural selection are a response 
to selective pressures that are often unknown. Estimating the age of selection signals may allow reconstruction 
of the history of environmental changes that shaped human phenotypes with specific ages associated with 
human dispersal out-of-Africa and the spread of agriculture (Nakagome et al., 2016). 
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Both mtDNA and Y-chromosomes have been used to reconstruct human history. The former (mtDNA) is 
often the better option for analyzing ancient DNA because it is easier to obtain, is present in higher number 
in human cells, and does not undergo recombination. However, mtDNA reflects only the maternal history of 
a population and the history of a single individual may not accurately reflect the history of a population. For 
this reason, mtDNA studies should be complemented by data on the male-specific Y chromosome 
(Pakendorf & Stoneking, 2005). Rates of human evolution have been estimated by applying a number of sub- 
stitution models to mtDNA or Y-chromosome sequence data (Behar et al., 2008; Kivisild, 2015; Soares et al., 
2013; Wang et al., 2014). In some examples, calibration has relied on the assumption that the genetic separa- 
tion between humans and chimpanzees occurred at ~6 Ma and that the evolutionary process has been clock- 
like since that time, while more complex models involve different substitution rates for coding and control 
regions of the genome (see Endicott et al., 2009). Molecular clocks can be tested by comparison with archae- 
ological data pertaining to human migration. 


Phylogeny based on the first complete mtDNA sequence data available in year 2000 (from 52 modern indi- 
viduals selected from around the world) yielded a time to the most recent common ancestor (TMRCA) of 
~190 ka (seeOppenheimer, 2003, 2009 ; Soares et al., 2009), apparently consistent with the conclusions of 
Cann et al. (1987) and with dated fossil finds of early AMH in Ethiopia (McDougall et al., 2005). More recent 
studies of human evolution from mtDNA (Oppenheimer, 2012), and of paternal evolution from Y chromo- 
somes (Poznik et al., 2016), yielded a broadly consistent picture of human evolution over the last ~200 kyr 
although the timing of branches in the evolutionary tree is poorly constrained and depends on estimated 
rates of mutation and population size (Figure 4). As pointed out by Wang et al. (2014), Y-chromosomal sub- 
stitution rates obtained using different calibration modes vary considerably and produce disparate recon- 
structions of human history. An additional determinant of substitution rate is the efficacy of purifying 
selection, which in turn depends not only on the particular constraints of each chromosome but also on 
the long-term effective population size for each chromosome type (Elhaik et al., 2014). Furthermore, muta- 
tion rates may have changed during hominoid evolution (Scally, 2016). 


The TMRCA is not a unique number, but rather a probability distribution based on two fundamental 
assumptions: the number of mutations and the mutation rate. Estimates of TMRCA depend strongly on 
the substitution rate, and different results can be obtained by using different rates. For example, Mendez 
et al. (2013) estimated a very early date of 338 ka for the TMRCA of the Y-chromosome tree (A00) from a 
population of African-Americans. The authors explained this early age by either long-standing population 
structure among modern human populations or archaic introgression from unknown species into the ances- 
tors of modern humans in western Central Africa. However, other researchers (e.g., Wang et al., 2014) have 
pointed out that this ancient TMRCA can be partially attributed to the low substitution rate used by the 
authors. By using either a higher mutation rate or more extensive sequencing data, the estimate of 
TMRCA becomes much younger at about 208 ka (Elhaik et al., 2014) or 257 ka (D'Atanasio et al., 2018). 


Poznik et al. (2013) reported entire Y-chromosome and mitochondrial genome sequences using a within- 
human calibration point to estimate the substitution rate. The results indicate small TMRCA differences 
from Y chromosomes (120-156 ka) and mtDNA (99-148 ka) that disagree with the conventional suggestion 
that the common ancestor of male lineages lived considerably more recently than that of female lineages. 
Analyzing the entire Y-chromosome dataset of the 1,000 Genomes Project (1,000 Genomes Project et al., 
2012; 1,000 Genomes Project, 2015) using the pedigree-based substitution rate, Wang et al. (2014) estimated 
the TMRCA at 105 ka that is consistent with the estimate (105 ka) of Cruciani et al. (2011) and the estimate 
(101-115 ka) of Wei et al. (2013). 


Phylogeny based on mtDNA sequence data has provided an estimate of TMRCA (190 ka) with branching 
episodes at ~120, ~70, ~40, and ~15 ka, which coincide with AMH dispersal patterns (Oppenheimer, 
2012). Wei et al. (2013) used five substitution models to assess phylogenetic nodes from the TMRCA of 29 
Y chromosomes that yielded branching ages for one of the models (GENETREE-2) of 112 +12 ka, 68 +7 
ka, 49 +6 ka, and 13 +2 ka. Applying the BEAST method to 68 worldwide Y chromosomes, Scozzari et al. 
(2014) dated the first two splits in their tree at ~196 and ~167 ka, followed by TMRCAs at ~110 ka, 85-77 
ka, 51-33 ka, and 8-22 ka (Figure 4). 


It remains challenging to reconstruct population structure prior to ~60 ka using existing data, but the 
modern-human African mtDNA pool contains phylogenetic patterns that can be used to estimate the ages 
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of several haplogroups. Analyzing LO HVS-I sequences, Rito et al. (2013) classified five branches (L0a, LOb, 
LOd, LOf, and LOk) of the mtDNA tree. The age estimates indicate that the mtDNA tree split to form LO at 
~180 ka, and later diversity follows a geographical distribution from southern Africa northward. Around 
~128 ka, two distinct AMH groups coexisted in Africa with a first split around 119 ka (LOk) followed by major 
clades at 98.7 ka (LOf), 70.9 ka (LOb), and 42.4 ka (L0a). It is noteworthy that the splitting of the widespread 
and common L0a, the diversification of LOk, and population increases of haplogroups LO and L3 can be 
dated to ~40 ka (Figure 4). 


8. Discussion 


Past episodes of low geomagnetic dipole field strength affect UVR levels arriving at the Earth's surface 
because reduction in magnetospheric shielding results in lower stratospheric ozone levels, and hence low- 
ered UVR shielding (Wei et al., 2014). Solar storms enhance nitrogen oxide production in the stratosphere 
that, in early 2004, led to ozone reduction of more than 60% at high northern latitudes (Randall et al., 
2005). Similarly, atmospheric NO concentrations are enhanced during times of low geomagnetic field 
strength when shielding from solar storms and GCRs is diminished. Earth's surface naturally emits N20 
from the oceans and from soils, with emissions having increased due to anthropogenic practices. N20 emis- 
sions are enhanced during interglacial and interstadial (warm) climate states, and increased by ~50% at the 
last glacial termination (Schilt et al., 2013, 2014). Enhanced atmospheric N2O concentrations during the 
Belling-Allerod warm period, which coincided with an apparent geomagnetic field strength minimum at 
~13 ka, would have elevated UVR reaching the Earth's surface at this time. It is noteworthy that the drastic 
magnetic field intensity decrease after the peak at 18 ka (Figure 1) coincided with paleoclimatic changes 
especially in the southern hemisphere (Boex et al., 2013; Moreno et al., 2015; Martinez-Garcia et al., 
2014). Evidence from Antarctic ice cores indicated sudden enhanced tropospheric UVR that has been related 
to a Mount Tahahe eruption dated to 17.7 ka (McConnell et al., 2017). 


A selection process resulted in Neanderthal disappearance at ~41 ka that apparently did not affect AMHs. 
Previous hypotheses for Neanderthal disappearance, and expansion of AMHs, include differences in subsis- 
tence strategies, language skills, and technical, economic, and social systems, and the ability to adapt to 
changing environments (e.g., Kochiyama et al., 2018). Response to environmental change such as UVR flux 
at Earth's surface would have involved the AhR, a chemosensor that regulates immunity and differs in 
AMHs versus Neanderthals and other primates (Hubbard et al., 2016). Ages for the end of the Mousterian 
tool industry and Neanderthal demise (Higham et al., 2014) are now tightly constrained to the Laschamp 
magnetic excursion (at 41 ka) implying a role for high UVR levels during the Laschamp field intensity 
minimum (Figure 3). 


Prominent low geomagnetic field intensity episodes at 285 ka, 190 ka, 100-120 ka, 64 ka, 41 ka, and 13 ka 
(Figure 1) appear to correspond to important times in evolution of hominins and other large mammals. 
At the LQE, megafauna in Australia, Europe, and the Americas were thought to have become extinct over 
a protracted time in the Late Quaternary; however, improved age estimates and discovery of new fossil sites 
(Barnosky & Lindsey, 2010; Cooper et al., 2015; Faith, 2014; Faith & Surovell, 2009; Koch & Barnosky, 2006; 
Price et al., 2011; Roberts et al., 2001; Saltré et al., 2016; Stuart, 2015; Wroe et al., 2013) have led to extinction 
peaks becoming progressively constrained to ~13 ka in the Americas, to ~40 ka (with earlier episodes at 84, 
107, and 122 ka) in Australia, and a complex combination of Late Pleistocene ages (including ~13 and ~40 
ka) in Eurasia and Africa. 


One of the outstanding and intriguing aspects of the LQE is the strong correlation between extinction and 
body mass. As the vast majority of cell mutations are deleterious, large long-lived organisms are at an evolu- 
tionary disadvantage. In addition, small mammals often have opportunities to avoid UVR through burrow- 
ing. In modern mammal populations, there is, however, no apparent correlation between body mass and 
cancer occurrence, known as Peto's Paradox (Peto et al., 1975). The elephant genome includes 20 copies of 
an oncogene (TP53) that is a crucial tumor suppressor gene involved in apoptosis in response to DNA 
damage, whereas other mammals usually have small numbers of this gene (Abegglen et al., 2015; Sulak 
et al., 2016). This discovery may explain why elephants are one of few large mammals to pass the LQE barrier. 
On the other hand, the remains of two extinct mammoth species included more than a dozen copies of TP53 
in their genomes, and the American mastodon, that disappeared at ~ 50 ka, had 3-8 copies in its genome 
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which implies an evolutionary selective trend in the increased number of copies of TP53 in the genomes of 
these related megafauna (Sulak et al., 2016). TP53 is considered the guardian of the genome due to its role 
in mitigating DNA damage, and is itself a target of UV-induced mutations (Aylon & Oren, 2011; de Pedro 
et al., 2018). 


Huang et al. (2017) analyzed the body weight of species in two orders of large ungulate herbivores 
(Artiodactyla and Perissodactyla) from the Neogene fossil record in Europe and North America. They found 
a significant and progressive increase in body weight from early Miocene to late Pliocene for both orders in 
North America and for Artiodactyla in Europe. This work was followed by the analysis of Smith et al. (2018) 
who documented a global increase in body-size of megafauna during the Cenozoic, with abrupt downsizing 
at the LQE in Late Pleistocene. We speculate that large mammals may have reached a natural body-size limit 
by Late Pleistocene, due to increased likelihood of cell mutation at times of high UVR flux, as life spans and 
body weights increased during the Cenozoic. The envisaged role of UVB in reducing megafaunal popula- 
tions at the LQE does not involve an instantaneous “blitzkrieg,” but rather an accumulation of UVB- 
triggered mutations over multiple (~30) generations or ~0.5-1 kyr, the approximate duration of RPI minima 
associated with magnetic excursions. 


There is no clear correlation of LQE events, either at ~13 ka (in North America and Eurasia) or ~40 ka (in 
Australia and Eurasia), with the first appearances of humans, which apparently preceded the LQE by at least 
10 kyr in Australia, and at least several kyr in North America. The human population of Australia at ~40 ka 
was likely no more than a few tens of thousands, with no evidence of an increase in population at this time 
(Williams, 2013). According to Webb (2013), the “overkill” hypothesis “is more sensational, granted, but the 
arguments are unrealistic and the evidence for it, at least in Australia, is non-existent.” The case for “over- 
kill” in North America at ~13 ka involves the close coincidence of the LQE with the brief (few century dura- 
tion) Clovis-tool horizon. On the other hand, based on Sporormiella decline in lake sediments from the 
eastern United States, the demise at North American herbivores began at ~14.5 ka (Gill et al., 2009, 2012), 
prior to the first appearance of Clovis tools. Clovis tools have only rarely been found in association with 
megafaunal remains. Evidence for pre-Clovis human occupation in the Americas includes locations in 
southern Chile (Dillehay et al., 2015), and the Florida panhandle where human occupation at ~14.5 ka pre- 
dated the Sporormiella decline by ~2 kyr (Halligan et al., 2016), consistent with Sporormiella disappearance 
at ~12.7 ka at another Florida panhandle location (Perrotti, 2018). Paucity of evidence for human occupation 
prior to ~13 ka can be attributed to habitation being concentrated in coastal regions that were largely sub- 
merged during the subsequent (last) deglaciation, when sea level rise necessitated migration of humans 
and other terrestrial mammals into the continental interior. We speculate that large mammals, particularly 
those that did not burrow, were particularly susceptible to DNA damage associated with low geomagnetic 
field strength and the resulting increase in UVB reaching the Earth's surface. The coincidence of low mag- 
netic field strength with LQE events, particularly at ~40 and ~13 ka, but also at 107 ka and 122 ka in 
Australia, implies that UVB flux was a contributing cause of the LQE in North America, Europe, 
and Australia. 


The ~190-ka paleointensity minimum at the Iceland Basin excursion corresponds to TMRCAs determined 
from the mtDNA and Y chromosomes of modern humans (Gonder et al., 2007; Oppenheimer, 2009, 2012; 
Poznik et al., 2016; Soares et al., 2009; Wei et al., 2013). Other paleointensity minima at 100-120 ka, 64 ka, 
41 ka, and 13 ka (Figure 1) correspond to branches in phylogeny estimated from mtDNA and Y-chromosome 
analyses (Figure 4). For example, Wei et al. (2013) used five models to assess phylogenetic nodes from the 
TMRCA of 29 Y chromosomes that yielded branching ages for one model (GENETREE-2) of 112 +12 ka, 
68 +7 ka, 49 +6 ka, and 13 +2 ka, corresponding closely to minima in RPI records (Figures 1 and 4). 
Although the transition from Middle Paleolithic to Upper Paleolithic at ~40 ka did not correspond to signif- 
icant changes in AMH anatomy in the fossil record, it is an important time for human population structure 
(Figure 4; Lahr, 2016) with mtDNA implying rapid population growth at a time of major advance in tool 
design and the first appearance of artwork (Stoneking, 1994). 


Phylogeography, the study of human dispersal, demonstrates that when migration occurs from one region to 
another, new mutations unique to that region accumulate (Soares et al., 2016). Local adaption to different 
habitats, including changes in exposure to mutagenic solar radiation, partially controlled by the magnetic 
field, are potential sources of phenotypical divergence (Jablonski & Chaplin, 2000, 2010). During 


CHANNELL AND VIGLIOTTI 


123 


ADVANCING EARTH 


Reviews of Geophysics 10.1029/2018RG000629 


migrations, humans have adapted to differences in climate, altitude, and resource availability. Migration to 
new environments alters selection pressures on the human genome, and genetic studies have identified cer- 
tain loci that were likely targets of selection for adaption to UVR (Ding et al., 2014). Highly pigmented skin 
protects against skin cancer but reduces synthesis of vitamin D3, so differences in the amount of UVR place 
differing selection pressures on pigmentation genes (Creanza & Feldman, 2016). 


9. Conclusions and Outlook 


Although coincidence does not prove causality, the timing of geomagnetic field strength minima (hence 
enhanced UVR flux at Earth's surface) appears to correspond to events in mammalian evolution. 
Improvements in knowledge of past geomagnetic field strength, new mammalian fossil finds, advances in 
radiocarbon dating and DNA analyses of fossils, use of dung-fungal proxies for herbivore populations, and 
advances in the use of mtDNA and Y chromosomes to map human phylogeny have all contributed to this 
possible linkage. Minima in geomagnetic dipole field strength at ~13 and ~41 ka (Laschamp; Figure 1) led 
to stratospheric ozone depletion and UVB levels at the Earth's surface that may have contributed to extinc- 
tion of large mammals at these times, although extinction dates are associated with large errors due to the 
inadequacy of the fossil record (Figure 2). The Neanderthals were apparently victims of the ~41 ka 
(Laschamp) magnetic field minimum, and differences in the AhR of modern humans and Neanderthals 
may explain why we passed the Laschamp evolutionary barrier. Episodes of low field strength further back 
in time notably at ~64 ka, ~100-120 ka (Blake; Figure 1), and ~190 ka (Iceland Basin; Figure 1) may have 
contributed to phylogenetic nodes in hominin evolution revealed by fossil finds and by studies of mtDNA 
and Y-chromosomes (Figure 4). According to this hypothesis, the geomagnetic field influenced evolution 
of large long-lived mammals through exposure to UVR at times of low field strength, with foci of the extinc- 
tion (e.g., Australia and Europe at ~41 ka, and North America and Europe at ~13 ka) depending on the spe- 
cific geometry of stratospheric ozone depletion during episodes of low field strength. Ozone holes are 
preferentially located at high latitudes because of the role of stratospheric temperatures and polar strato- 
spheric clouds in ozone depletion. UVR arriving at the Earth's surface may have had an influence on evolu- 
tion due its strong mutagenic effects, its potential for promoting oxidative damage on membranes and 
proteins, and the role of AhR in translocating UVB stress response to the nucleus. Lower levels of UVR 
reaching Earth's surface at low latitudes, due to the role of polar-stratospheric clouds and stratospheric tem- 
perature in ozone depletion, may partly explain the relative diversity of modern African megafauna. 


As the chronologies of both the RPI record and the paleontological record become better resolved, we spec- 
ulate that unconvincing efforts to establish a linkage between polarity reversal and extinction over the last 
~50 years will, in future, be superseded by a linkage between extinction and paleointensity minima particu- 
larly for land-dwelling mammals at higher latitudes where the magnetic field strength is linked to loss of 
stratospheric ozone and enhanced UVR reaching Earth's surface. The importance of Australian extinction 
events at the time of Laschamp excursion (41 ka) and the Blake excursions (100 and 120 ka) may imply that 
stratospheric ozone depletion was largely in the Southern Hemisphere at these times, and in the Northern 
Hemisphere when the North American extinctions are manifested at ~13 ka. 


The magnetic dipole moment has decreased by ~10% since 1833 (Gauss' first direct field intensity measure- 
ment) or ~5% per century. Although the present dipole field strength may not be appreciably lower than the 
average during the Brunhes Chron, the current field intensity decrease, combined with asymmetry in the 
field detected in satellite data, has led to speculations that the geomagnetic field may reach intensity levels 
appropriate for a magnetic excursion or polarity reversal in the next 1,000-2,000 years (Glassmeier et al., 
2009; Hulot et al., 2002; Laj & Kissel, 2015). 


Twenty years ago, no more than five Quaternary magnetic excursions had been recorded (see Opdyke & 
Channell, 1996), whereas today ~20 Quaternary magnetic excursions have credible documentation (Laj & 
Channell, 2015). In the next few decades, we speculate that integration of excursion records with RPI and 
cosmogenic isotope records, documentation of mammalian extinction, understanding the role of geomag- 
netic field in UVR flux, and the role of Ahr in modulating the deleterious effects of UVR in mammals, as well 
as improvements in human phylogeny from mtDNA and Y chromosomes, will clarify the role of the geomag- 
netic field in mammalian evolution. 
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Glossary 


6-4 photoproducts 
(6-4 PPs) 


Anhysteretic 
remanent 
magnetization 
(ARM) 


Apoptosis 


Artiodactyla 


Aryl hydrocarbon 
receptor (AhR) 


Brunhes Chron 


Campanian 


Ignimbrite (CI) 


Chlorofluorocarbons 
(CFCs) 


Clades 


Clovis horizon 


Coercivity 


6-4 pyrimidine-pyrimidone photoproducts (6-4 PPs) and cyclobutane pyrimidine 
dimers (CPDs, see below) are the common ultraviolet radiation (UVR) products 
formed in the human skin during exposure to sunlight. Although the 6-4 PPs are 
not as prevalent as CPDs (in lesion formation), they are significantly more 
mutagenic. 

A laboratory-induced remanent magnetization imparted by a strong alternating 
magnetic field in the presence of a weak direct current bias field. The bias field 
usually has an intensity comparable to the geomagnetic field. ARM provides an 
indication of magnetic mineral concentrations, and the efficiency of ARM 
acquisition varies with magnetic grain size. 

A form of programmed cell death that occurs in multicellular organisms. It is a 
highly regulated and controlled process that confers advantages during an 
organism's lifecycle. Because apoptosis cannot be halted once initiated, it is a 
highly regulated process. 

Order of even-toed ungulates, hoofed animals that bear their weight on two toes 
including hippopotamuses, camels, antelopes, sheep, giraffes, and deer. 

A protein that regulates gene expression, cell physiology, and organ homeostasis 
with a key role in skin integrity and immunity. Its functions are related to cell 
proliferation, adhesion, and migration as well as to cell differentiation. It has a 
strong function in mediating the response to xenobiotic and toxins particularly 
Dioxins, and to UVR. 

The most recent time interval of the geomagnetic polarity timescale, 
characterized by a normal polarity of the geomagnetic field that has existed 
since ~773 ka. 

The CI eruption was located near Napoli (Italy) and was probably the largest 
volcanic eruption in the Mediterranean region during the last 200 kyr. CI ash 
was deposited on the Russian Plain, and throughout the Eastern Mediterranean 
and northern Africa. The event coincided with the onset of a cold climatic phase 
known as Heinrich Stadial 4 (HS4) approximately 40,000 years ago. 

An organic compound (commonly known as Freon) that contains carbon, chlorine, 
and fluorine, produced as a volatile derivative of methane and ethane. Many CFCs 
have been commonly used as refrigerants, propellants (in aerosol applications), and 
solvents. Although the concentration of CFCs in the atmosphere is very small, 
measured in parts per trillion, they do contribute significantly to the 
enhancement of the greenhouse effect. CFCs contribute significantly to ozone 
depletion in the upper atmosphere. The manufacture of such compounds has 
been phased out under the Montreal Protocol, and their use is being replaced by 
other products such as hydrofluorocarbons. 

In taxonomy, a clade is defined as a group of organisms consisting of a single 
common ancestor and all descendants of that ancestor. 

The short-lived Clovis horizon comprises distinctive stone tools marking 
early human occupation of North America at ~13 ka. The Clovis culture is a 
prehistoric Paleo-Indian culture, named after characteristic stone tools found in 
close association with Late Pleistocene fauna at Blackwater near Clovis, New 
Mexico, in the 1920s and 1930s. Most of the indigenous cultures of the Americas 
are considered to descend from the Clovis people. 

The coercivity, also called magnetic coercivity, is a measure of the ability of a remanent 
magnetization to withstand an external demagnetizing magnetic field. For 
ferromagnetic materials (sensu lato), the coercivity (or coercive force) is the intensity 
of the applied magnetic field required to reduce the remanent magnetization of that 
material to zero after the magnetization of the sample has been driven to saturation 
in the opposite direction. Thus, coercivity measures the resistance of a ferromagnetic 
material (sensu lato) to demagnetization by an external field. 
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Cosmogenic isotopes 


Cyclobutane 
pyrimidine 
dimers (CPDs) 
Cytosol 


Deoxyribosyl 


DNA photoproducts 


Endogenous ligand 


Energetic particle 
precipitation (EPP) 


Gaussian-resampled, 
inverse-weighted 
McInerny 

et al. (GRIWM) 


Geomagnetic 
excursions 


Radioactive isotopes created in the upper atmosphere when galactic cosmic rays 
collide with atmospheric molecules at high speed. The production rate of these iso- 
topes depends on the intensity of the cosmic radiation, which is related to the 
strength of the Earth's magnetic field and solar activity. Therefore, records of cos- 
mogenic isotope production are useful as a proxy for the strength of Earth's mag- 
netic field and for understanding the relation between the Earth magnetic field, 
and variations in solar activity. 

One of two major photoproducts generated by UV irradiation. CPDs are molecular 
lesions formed in the human skin during exposure to sunlight due to UV radiation 
and photochemical reactions. 

It is the fluid inside living cells. Proteins, organelles, and other structures of cells 
live in this water-based fluid. 

A univalent radical derived from deoxyribose, a monosaccharide (simple sugar). It 
is derived from the sugar as indicated by its name: deoxy sugar. Deoxyribose is the 
five-carbon sugar molecule involved to form the phosphate backbone of DNA 
molecules. DNA, or deoxyribonucleic acid, is a polymer formed of several nucleic 
acids. Each nucleic acid is composed by a molecule of deoxyribose bound to a 
phosphate group and either a purine or a pyrimidine. As part of DNA, 2- 
deoxyribose derivatives have a significant role in biology. The DNA molecule, 
which is the main repository of genetic information in life, consists of a long 
chain of deoxyribose-containing units called nucleotides, linked via phosphate 
groups. 

Exposure to the ultraviolet component of sunlight causes DNA damage, which 
subsequently leads to mutations, cellular transformation, and cell death through 
the creation of different photoproducts. There are two important DNA 
photoproducts, namely CPDs and pyrimidine-pyrimidone photoproducts (6-4 PPs). 
A ligand is a protein that attaches (binds) to another protein called a receptor; the 
latter have specific sites into which the ligands fit like keys into locks. Ligands that 
are produced in the body are called endogenous. 

Precipitation of a group of highly energetic electrons, protons, neutrons, and ions 
that are accelerated into the atmosphere through various heliophysical and 
geomagnetic processes. They enter the atmosphere mainly in the vicinity of the 
geomagnetic poles. When energetic particles enter the atmosphere they ionize 
and dissociate atmospheric constituents, resulting in the formation of reactive 
nitrogen oxides. EPPs have been shown to contribute up to 10% of the 
stratospheric NO, budget and up to 40% of the polar stratospheric NO, budget. 
Once in the stratosphere, NO, produced by EPPs (EPP-NO,) interferes with 
catalytic cycles involving ozone (O3). Theoretically, changes in O, can also lead 
to changes in temperature and winds, which means that EPPs have the potential 
to impact climate. 

Dating method for estimating the probability of extinction by using an approach 
that weighs observations inversely according to their temporal distance from the 
last observation of a species’ confirmed occurrence. For dates with associated 
radiometric errors it is able to sample individual dates from an underlying 
fossilization probability distribution. 

Excursions represent a significant, but short-lived, change of direction of the 
Earth's magnetic field and are apparently manifested globally. Intervals of 
constant polarity are punctuated by geomagnetic excursions where magnetic 
directions depart from the usual geocentric axial dipole, and when adequately 
recorded, achieve the opposite polarity direction for a short time. Excursions 
have been recorded in both volcanic and sedimentary records. These 
directional aberrations typically have durations of a few thousand years, or less 
than 1 kyr in some cases, and are characterized by a decrease in strength of 
the main axial dipole. 
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Halons 


Haplogroup 


Haplotype 


Heinrich Stadial 
(HS) 4 


Isothermal remanent 


magnetization (IRM) 


Late Quaternary 
extinction (LQE) 


Ligand binding 


Luminescence age 
dating 


Mediterranean 
sapropel S7 


A group of hydrocarbon compounds in which some hydrogen atoms are 
replaced by bromine and fluorine atoms, and sometimes also chlorine. Halons 
are nonconductors of electricity and may be utilized in fighting fires in ignitable 
liquids and most solid flammable materials. Halons in the atmosphere are 
responsible for ozone depletion and contribute to greenhouse warming. 

Group of individuals sharing a common ancestor identified by DNA sequences 
defined by shared mutations and which tend to show regional specificity. In 
human genetics, the most commonly studied haplogroups are Y chromosome (Y- 
DNA) haplogroups and mitochondrial DNA (mtDNA) haplogroups, each of 
which can be used to describe genetic populations. 

A group of genes in the chromosome of an organism that are inherited together, 
because of genetic linkage from a single parent. The word “haplotype” is derived 
from the word “haploid” that describes cells with only one set of chromosomes 
and from the word “genotype,” which refers to the genetic makeup of an 
organism. 

A Heinrich event is a natural phenomenon. Several occured during a period 
known as Marine Isotope Stage 3 (about 59-27 ka). Large icebergs 
containing rock mass eroded by the glaciers broke off from glaciers and 
traversed the North Atlantic, as they melted this material was deposited on 
the sea floor as ice rafted debris. One of the stronger Heinrich stadials is 
HS4 at ~40 ka. 

IRM is an artificial (usually laboratory acquired) remanent magnetization acquired 
by applying a strong direct field to a ferromagnetic (sensu lato) material. It is a 
useful parameter used for detection of magnetic minerals with high coercivity 
such as hematite or goethite. 

Widespread extinction of large terrestrial mammals (megafauna) during the Late 
Quaternary, observed in the Americas, Europe, Africa, Asia, and Australia. As 
the ages of extinction have become better constrained, the timing of these 
extinctions has become focused in the vicinity of 40 and 13 ka. 

Intermolecular interactions occurring between or among proteins, nucleic acids, or 
small molecules. The interaction of ligands with their binding sites is often 
characterized in terms of a binding affinity. In general, high-affinity ligand 
binding occurs from greater intermolecular force between the ligand and its 
receptor, whereas low-affinity ligand binding involves less intermolecular force 
between the ligand and its receptor. 

Refers to a group of dating methods based on the determination of how long ago 
mineral grains were last exposed to sunlight. The techniques include optically 
stimulated luminescence, infrared stimulated luminescence, and 
thermoluminescence dating. It is useful to geologists and archaeologists who 
want to know the age of a particular event. Sediments and soils contain traces of 
radioactive isotopes such as potassium, uranium, thorium, and rubidium. The 
ionizing radiation emitted during the decay of these elements is absorbed and 
trapped by quartz and potassium feldspar grains occurring in the sediment. 
Stimulating these mineral grains using either light (blue or green for optically 
stimulated luminescence; infrared for infrared stimulated luminescence) or heat 
(for thermoluminescence dating) produces a luminescence signal emitted as 
stored unstable electron energy, the intensity of which varies depending on the 
amount of radiation absorbed during burial, and specific properties of the mineral. 
The sediments of the Mediterranean (central-eastern) are characterized by the 
presence of organic-rich (>2% organic carbon) layers called sapropels. They were 
deposited during periods of reduced oxygenation related to increasing riverine 
input occurring in correspondence with peaks in solar insolation. Each peak has 
been numbered starting from the youngest (S1; 8-10 ka) and the sapropel S7 
occurred around 190 ka. 


CHANNELL AND VIGLIOTTI 


727 


ADVANCING EARTH 


Reviews of Geophysics 


10.1029/2018RG000629 


Melanocyte homeos- 
tasis 


Mitochondria 


Mitochondrial DNA 
(mtDNA) 


Mitochondrial Eve 


Mitohormesis 


Mousterian tool 
industry 


Natural remanent 
magnetization 
(NRM) 


Pedigree-based 
substitution rate 


Perissodactyla 


Peto's paradox 


The maintenance of systems within a cell as conditions change is 
called homeostatic regulation, and melanocyte (a mature melanin-forming skin 
cell) homeostasis is a paradigm for understanding the formation of melanoma. 
Mitochondria are part of a eukaryote cell found in the cytoplasm. They 
oxidise glucose to provide energy for the cell through the creation of a molecule, 
called adenosine triphosphate that cells use as an energy source. 

Mitochondria are tiny organelles that live in the cytoplasm of cells, 
called mitochondria. Each cell contains thousands of mitochondria with its own 
small circle of DNA, a reminder of their distant bacterial ancestry. MtDNA is not 
located in the cell nucleus, and in most species, including humans, it is inherited 
solely from the mother and is not subject to genetic recombination during 
meiosis. Therefore, it remains unchanged from generation to generation. Human 
mtDNA was the first significant part of the human genome to be sequenced. 
Since animal mtDNA evolves faster than nuclear genetic markers, it represents a 
pillar of phylogenetics and evolutionary biology. It also makes possible an 
examination of the affinity of populations, and so has become important in 
anthropology and biogeography. 

Analyses of the mitochondrial DNA of living humans have shown that we descend 
from a common female ancestor that has been dated to about 200 ka. This woman 
has been named “Eve” or “Mitochondrial Eve” after a journalist's confusing 
reference to the unrelated biblical story of the first woman being created by God 
from Adam's rib. (“The man called his wife Eve [“life”], because she was the 
mother of all who live.” Genesis 3:20) 

Mitohormesis, known also as mitochondrial hormesis, is a particular form of 
hormesis that is a nonlinear response to potentially harmful substances. 
Mitochondrial function often results in excessive production of reactive oxygen 
species (ROS) that are responsible for many chronic diseases. However, moderate 
levels of mitochondrial ROS can protect against chronic disease by stimulating 
mitochondrial capacity and endogenous antioxidant defense. This phenomenon 
is called mitohormesis. 

The predominant industry of the Middle Paleolithic is termed Mousterian, named 
for its type location at Le Moustier, in Dordogne, France. Mousterian industry is the 
tool culture traditionally associated with Neanderthal man in Europe, western 
Asia, and northern Africa. Mousterian tools disappeared abruptly from Europe 
with the passing of Neanderthal man. 

The magnetic remanence carried by rocks or sediments before the laboratory 
demagnetization or magnetization treatments. The NRM of a rock or sediment is 
usually represented by multiple magnetization components acquired through its 
geologic history. Thermal or alternating field laboratory treatments are usually 
necessary to identify the primary remanence acquired at the time of rock 
formation or sediment deposition. 

A method to estimate mutation rate by comparing the mtDNA sequences of 
a sample of parent/offspring pairs or analyzing mtDNA sequences 
of individuals from a deep-rooted (well established) genealogy. The number 
of new mutations in the sample is counted and divided by the total 
number of parent-to-child DNA transmission events to calculate a mutation 
rate. 

Order of odd-toed ungulates, hoofed animals that bear most of their weight on a 
single toe including horses, zebras, rhinoceroses, and tapirs. 

Peto's paradox, named after Richard Peto, is the observation that at species level, 
the incidence of cancer does not appear to correlate with the number of cells in 
an organism (e.g., body weight). Although cell mutations are usually deleterious, 
the incidence of cancer is not apparently related to the body-size of individual 
mammals. 
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Phenotypical 


Phylogeny 


Polar-stratospheric 
clouds 


Quaternary 


Reactive oxygen 
species (ROS) 


Relative 
paleointensity 
(RPI) 


Signor-Lipps effect 


Sclerophyll 
vegetation 


Solar proton events 
(SPE) 


The observable physical or biochemical characteristics of an organism determined 
by genetic and environmental influences. 

The evolutionary history and relationships among people or groups of 
organisms (e.g., species or populations). These relationships can be discovered 
through phylogenetic methods that evaluate observed heritable traits, such as 
DNA sequences or morphology, using an evolutionary model for these traits. 
The results of these analyses represent the phylogeny (also referred to as the 
phylogenetic tree). 

Polar-stratospheric clouds are typically formed in the stratosphere (15- to 25-km 
altitude) at high latitudes especially in winter because they require very low 
temperatures. They are implicated in the formation of ozone holes because they 
convert benign forms of chlorine into reactive, ozone-destroying forms, and they 
remove gaseous nitric acid that would otherwise moderate the destructive impact 
of chlorine. 

The most recent of the three periods of the Cenozoic Era in the geologic time scale. 
It spans the last 2.588 Myr and is divided into two epochs: the Pleistocene (2.588 
million years ago to 11.7 thousand years ago) and the Holocene (11.7 thousand 
years ago to present). 

Oxygen-derived molecules that act as powerful oxidants and may form via a 
large number of physiologic and nonphysiologic processes as a result of 
natural consequences of aerobic metabolism. Although these molecules play 
a role in the oxygen-dependent defense mechanism against bacteria, they 
may also be highly damaging, as they can attack biologic macromolecules, 
including lipids, proteins, and DNA, and lead to significant tissue damage. 
In normal conditions, ROS can act as immune system modulation and can 
activate various signal transduction pathways. ROS levels can increase 
dramatically during times of environmental stress (e.g, UVR or heat 
exposure). 

The intensity of the Earth's axial dipole field is globally coherent and can be 
reconstructed in the past by studying suitable rock and sediment samples or 
archaeological materials. The RPI proxy in sediments usually comprises the NRM 
intensity normalized by the intensity of a laboratory-acquired magnetization 
designed to activate the same population of magnetic grains that carry the NRM, 
thus compensating for variations in concentration of NRM-carrying grains 
through the sedimentary section. The laboratory-applied normalizers are 
typically ARM and/or IRM. 

The first and last occurrence of a taxon will never be accurately recorded in the 
fossil record because of the inherent incompleteness of the record. 

Type of vegetation characterized by hard, leathery, evergreen foliage well 
adapted to prevent moisture loss. Broad-leaved sclerophyll vegetation, 
including species such as holly (Ilex), is known as Mediterranean vegetation 
because it develops in regions with a Mediterranean climate—hot, dry 
summers and mild, wet winters. Pines are examples of narrow-leaved 
sclerophyll vegetation. 

ASPE or “proton storm” occurs when particles (mostly protons) emitted by the Sun 
are accelerated either close to the Sun during a flare or in interplanetary space by 
coronal mass ejection. Solar protons normally have insufficient energy to 
penetrate the Earth's protective magnetosphere. However, during unusually 
strong flares, protons can be produced with sufficient energies to penetrate the 
Earth's magnetosphere and ionosphere around the poles. They can pass through 
the Earth's magnetic field and cause ionization in the ionosphere. The effect is 
analogous to an auroral event, except that protons instead of electrons are 
involved. Energetic protons represent a significant radiation hazard to spacecraft 
and astronauts. 
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The time to the most Time to the most recent individual from which an entire population of the group is 
recent common directly descended. The term is used in reference to the ancestry of groups of genes 
ancestor (TMRCA) (haplotypes) or species. 

TP53 (p53) gene Tumor protein 53 is a protein-coding gene that acts as a tumor suppressor for many 
tumor types and induces growth arrest or apoptosis depending on physiological 
circumstances and cell type. The TP53 gene is responsible for instructions to 
create a protein called tumor protein p53. This protein acts as a tumor 
suppressor, which means that it regulates cell division by keeping cells from 
growing and dividing (proliferating) too fast or in an uncontrolled way. The p53 
protein is found in the nucleus of cells throughout the body, where it attaches 
(binds) directly to DNA. If the DNA in a cell is damaged by agents such as toxic 
chemicals, radiation, or UVR, this protein plays a critical role as to whether the 
DNA can be repaired or the damaged cell will self-destruct (undergo apoptosis). 
If it is possible to repair the DNA, p53 activates other genes to fix the damage. On 
the other hand, this protein can prevent the cell from dividing and signals it to 
undergo apoptosis. By stopping cells with mutated or damaged DNA from 
dividing, p53 helps prevent tumor development. 

Tryptophan Tryptophan is a a-amino acid that is used in the biosynthesis of proteins. Many 
animals (including humans) cannot synthesize tryptophan: they obtain it 
through their diet. Tryptophan is among the less common amino acids found in 
proteins, but it has critical structural and functional roles. 

Virtual axial dipole Intensity of an imaginary axial (along the Earth's rotation axis) centric (located in 
moment (VADM) the center of the Earth) dipole that would produce the estimated magnetic field 
reconstructed from paleomagnetic data. The reconstruction is virtual because the 

VADM is an approximation of a geocentric axial dipole moment. 

Xenobiotics A xenobiotic is a chemical substance within an organism that is foreign to the 
biological system and not naturally produced or expected to be present within the 
organism. It can also cover substances that are present in much higher 
concentrations than usual. The term is often used in the context of pollutants and 
their effect on biota. 

Y chromosome One of the two sex chromosomes in humans (the other is the X chromosome) and 
most mammals. The Y chromosome occurs in males, who have one X and one Y 
chromosome, while females have two X chromosomes. The Y chromosome 
contains genes that provide instructions for making proteins. Because only males 
have the Y chromosome, the genes in this chromosome are involved in male sex 
determination and development. The Y chromosome represents almost 2% of the 
total DNA in male cells. 
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Erratum 


In the originally published version of this article, discoloration appeared in Figure 2 and Figure 4 was not 
large enough. Figures 2 and 4 have been corrected, and this version may be considered the authoritative 
version of record. 
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